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Within the QCD factorization approach, we calculate the branching fractions and CP asymme- 
try parameters of 12 B + Kò (1430) Kk (+) decay modes under the assumption that the scalar meson 
kK (1430) is the first excited state or the lowest lying ground state in the quark model. We find 
that the decay modes with the scalar meson emitted, have large branching fractions due to the 


enhancement of large chiral factor ry °. The branching fractions of decays with the vector meson 


emitted, become much smaller owing to the smaller factor tr i 


. Moreover, the annihilation type 


diagram will induce large uncertainties because of the extra free parameters dealing with the end- 
point singularity. For the pure annihilation type decays, our predictions are smaller than that from 
PQCD approach by 2-3 orders of magnitudes. These results will be tested by the ongoing LHCb 
experiment, forthcoming Belle-II experiment and the proposing circular electron-positron collider. 


I. INTRODUCTION 


Although the quark model has made great success in 
describing most of the hadronic states, the lowest lying 
scalar mesons are too light to fit in the quark model. 
Two possible scenarios have been proposed on the ba- 
sis of whether the scalars lower than 1 GeV belong to 
the four-quark states or the classical two quark states. 
They are controversial for decades [1]. In scenario-1 
(S1), the scalars such as «(800), ao(980) and fo(980) are 
naively seen as the lowest lying qq states, and Kj (1430), 
ao(1450), and fo(1500) are the first excited states, cor- 
respondingly. In contrast, Aj (1430), ao(1450), (1500) 
are treated as the the qq ground states and their first 
excited states are about (2.0 ~ 2.3) GeV in scenario-2 
(S2), in which the lightest scalar mesons are the four- 
quark bound states. 


In the past decade, many B decay modes involving 
scalar mesons have been reported by BaBar, Belle and 
large hadron collider-b (LHCb) experiments. This may 
provide a unique feature to distinguish these two scenar- 
ios by the B meson tag [2]. It is thus hoped that the 
combination of the precise experimental measurements 
and the accurate theoretical predictions might provide 
us valuable information on the nature of scalar mesons. 


To achieve this goal, some hadronic B,(q = u, d, s) de- 
cays to scalar mesons have been studied in detail in the 
framework of the QCD factorization (QCDF) [1, 3, 4] 
and the perturbative QCD approach (PQCD) [5, 6]. Us- 
ing the QCDF approach, H-Y Cheng et.al had calculated 
the branching fractions and direct CP asymmetries of 
most decay modes in refs. [3, 4], such as B > fok, 
B —> K(1430)¢(p) and B > Kğ(1430)r, where most 
results accommodate the data with large uncertainties. 
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Very recently, the decays B + K§(1430)K) dominated 
by b > d penguin operators, have been calculated within 
the PQCD approach [6], however they have not been 
touched in QCDF frame. To complete, we therefore shall 
calculate the branching fractions of B + K§ (1430) Kk) 
decays in QCDF, as well as their CP asymmetry param- 
eters. In the experimental side, these observables are 
too small to be measured at current experiments, but 
they are hopeful to be measured in the future experi- 
ments such as the updated LHCb experiment, the high 
luminosity Belle-II experiment and the proposing high 
energy circular electron-positron collider. 

The paper is organized as follows. In Section 2, we 
give the analytic formula including the effective Hamil- 
tonian, the form factor and all corrections to the ampli- 
tudes. Presentation of results and discussions are given 
in Section 3. At last, we summarize this work in Section 
4. 


II. ANALYTIC FORMULA 


In this section, we shall start from the weak effective 
Hamiltonian responsible for B > K}(1430)K“ decays. 
In the standard model, it could be written as [7] 
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The explicit form of the operators O; and the corre- 
sponding Wilson Coefficients C; at different scale u could 
be found in ref.[7]. Vice» and Vaja are the Cabibbo- 
Kabayashi-Maskawa (CKM) matrix elements. Note that 
O;,2 are tree operators and others O3—10,7y,8g are pen- 
guin ones. 


In dealing with the nonleptonic charmless B decays, 
the decay amplitude is usually separated into the emis- 
sion part and the annihilation part in terms of the struc- 
ture of the topological diagrams. According to the fac- 
torization approximation based on the heavy quark limit, 
the former part could be written as the product of decay 
constant and form factor. While for the latter one, it is 
always regarded as being power suppressed. In QCDF [8] 
based on collinear factorization, the contribution of the 
non-perturbative sector is dominated by the form factors, 
and the non-factorizable effect in the hadronic matrix el- 
ements is controlled by hard gluon exchange. Thus, the 
hadronic matrix elements of the decay can be written as 
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where T, H and T/! denote the perturbative short- 
distance interactions, which can be calculated perturba- 
tively. x(x) (X = B, Mı, M2) are the universal non- 
perturbative light-cone distribution amplitudes that can 
be estimated by the light-cone QCD sum rules. The form 
factors of B meson to the pseudoscalar (P), the vector 
(V) and scalar (S) mesons transitions F?~™: being a 
part of the nonperturbative sector of QCD, lack a pre- 
cise solution. For the B > P and B —> V transition 
form factors, we will employ the results of the QCD sum 
rule method [9], since they have been used in calculating 
B = PP,PV and VV modes widely. As the form fac- 
tor of F?>*% as concerned, to the best of our knowledge, 
a number of approaches had been advocated to calcu- 
late them, such as QCD sum rule [10], light-cone QCD 
sum rule [11], PQCD [12] and covariant light front quark 
model (cLFQM) [13]. In this work, we use the results 
obtained by cLFQM [13] for keeping consistent with the 
results of other B > SP(V) decay modes [3]. 

Following the standard procedure of QCDF approach, 
the emission part of decay amplitude could be written as 


As(B° + M,M2) 
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(M,M2|O;|B) r is the factorizable matrix element, which 
can be factorized into a form factor times a decay con- 
stant, as stated before. The effective parameters a? can 
be calculated perturbatively, whose expressions are given 
by 
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with 7 = 1,--- ,10. The upper (lower) signs apply when 
i is odd (even), Cr = (N? — 1)/(2N.) with Ne = 3. 
The quantities V;(M2) account for vertex corrections, 
H;(Mı Mz) for hard spectator interactions with a hard 
gluon exchange between the emitted meson and the spec- 
tator quark of the B meson and P;(M2) for penguin con- 
tractions. Similarly, the annihilation contributions are 
described by the terms b;, and b; mw, which have the 
expressions 
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where the subscripts 1,2,3 of Af stand for the an- 
nihilation amplitudes induced from (V — A)(V — A), 
(V — A)(V + A) and (S — P)(S + P) operators, respec- 
tively. The superscripts i and f refer to gluon emission 
from the initial and final-state quarks, respectively. It 
should be stressed that the decays Bs K+K- and 
B >K ()+ KT are only induced by the annihilations 
type diagrams. Hereafter, the order of the a? (MM2) 
coefficients is dictated by the subscript Mı Mə, where 
Mə is the emitted meson and Mı shares the same spec- 
tator quark with the B meson. For the b;(MıMə2) of 
the annihilation, Mı means the one containing an anti- 
quark from the weak vertex, while Mə contains a quark 
from the weak vertex. The explicit expressions of Vj, 
H;,P;, bi, bi œw and their inner functions could be found 
in refs.[1, 3]. 

When dealing with the hard-scattering spectator and 
the weak annihilation contributions, we suffer from in- 
frared endpoint singularities X = i, dx/(1 — x), that 
cannot be calculated from the first principle in the QCDF 
approach and only be estimated phenomenologically with 
a large uncertainty. Following the arguments of ref.[8], 
we also parameterize these kinds of contributions by the 
complex quantities, Xy and X4 namely, 


Xy = (1+ one) n ZE, (6) 
Xa = (1+ paeta) n= , (7) 


where A = 0.5GeV. pu, pa are real parameters, and ¢# 
and ¢, are free strong phases in the range [—180°, 180°]. 
All the above four parameters should be fixed by the 
experimental data, such as branching fractions and CP 


asymmetries. In the so-called PQCD approach [14], one Including the emission and annihilation contributions, 
can eliminate these end-point singularities by keeping all the decay amplitude can be finally given as 
small transverse momenta of gluons and inner quarks. 
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II. NUMERICAL RESULTS AND DISCUSSION 


To proceed, we shall present numerical results obtained 
from the formulas given in the previous section. Firstly, 
we introduce the adopted parameters in our calculation. 
Secondly, we give the numerical results and show the the- 
oretical errors due to uncertainty of some parameters. At 
last, some discussions and comparisons will be added. 

The scalar meson, unlike the pseudoscalar one, has two 
kinds of decay constants, namely the vector decay con- 
stant fs and the scale-dependent scalar decay constant 
fs, the definitions of which are give by: 


(Ki (@)lEyud]0) = frs Pu 
(K9(p)|54)0) = mx Fr; (21) 
The two decay constants satisfy 


where m,(u) and malu) are the running current quark 
masses. It should be stressed that the decay constants 


m 


of Kõ(1430) have the signs flipped from S1 to S2. The 
definition of the form factors for the B —> S transitions 
are given by [13] 
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where P, = (p + p’), and q, = (p—p'),. In cLFQM , 
the momentum dependence of the form factor could be 
parameterized in a di-pole model form [13], 
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Together with the decay constants, the parameters F'(0), 
a and b for B —> S transitions in the different scenarios 
are summarized in Table I. 

The twist-2 light-cone distribution amplitude (LCDA) 
s(x) and twist-3 P(x) and PZ(x) for the scalar meson 


Ko made of the quarks sd are given by 
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TABLE I: The Parameters of the A (1430) in different scenarios 


Fo 7") a b |F 0) a b | Fixx (MeV) Bı Bs 


S1 0.26 
S2 0.21 


; 0.44 0.05 0.26 1.52 0.64|—300 +30 0.39 + 0.05 
; 0.44 0.05 0.21 1.52 0.64| 445 +50 —0.39 +0.09 —0.25+ 0.13 


—0.70 + 0.05 


TABLE II: Summary of input parameters 
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0.225 0.818 0.376 68° 250MeV 1.54ps 1.67ps 0.35 1/132 
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236MeV 5.36GeV 131MeV 221MeV 175MeV 0.49GeV 0.89GeV 1.43GeV 0.35 0.34 


with z = z2 — 21, f = 1 — x, and their normalizations are 
1 
f avbs(0) = fs, (26) 
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The above definitions of LCDAs can be combined into a 
single matrix element as 
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In general, the twist-2 light-cone distribution amplitude 
s has the form 


Ðs(z, u) = fs(u) 6z(1 — z) 
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where Bm(u) are scale-dependence Gegenbauer moments 
and Cy : (u) are the Gegenbauer polynomials. The Bm 
of different scenarios are also presented in Table I. As for 
the twist-3 distribution amplitudes, we shall adopt the 
asymptotic form for simplicity, shown as 


© (x) = fs, 


For the LCDAs of K‘*), we will employ the formulaes 
obtained from the QCD sum rules [15, 16]. The other 


Z(z) = fs 6z(1— z). (30) 
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used parameters, such as the CKM elements, the decay 
constants of the pseudo-scalar and the vector, and the 
form factors of B > K®), are also listed in Table II for 
convenience. 


Now, we turn to discuss the numerical results of the 
concerned decay modes. As we had stated in previous 
section, when calculating the hard spectator and the an- 
nihilation contributions, two endpoint singularities are 
parameterized phenomenologically by four free parame- 
ters, namely pH, QH and pa,pa, which should be deter- 
mined from the experimental data. In ref.[3], a global fit 
of pa and ġa to the B > SP data indicates pa = 0.15 
and ġa = 82° with x? = 8.3, so in this work for the 
central values (or “default” results), pa, = 0.15 and 
a,u = 82° are adopted. In Table III, we list the cal- 
culated branching fractions of B + Kë(1430)K®) de- 
cays, where the first uncertainty is due to the variations 
of Bı and fg, the second comes from the form factors 
and the strange quark mass, and the last one is induced 
by weak annihilation and hard spectator interactions in 
ranges pa. €E [0,0.3] and a,x € [0, 180°]. 


Here, we shall take B~ —> K- K§°(1430) and B7 > 
Ko (1430) K° as examples to illustrate the relative size of 
each contribution. The decay amplitude formula of each 
mode has been presented in eqs.(8) and (9), respectively. 
The first decay mode is characterized by B —> P transi- 
tion with scalar meson emitted; while the second decay 
mode is characterized by B > S transition with pseu- 
doscalar meson emitted. Because of the small vector de- 
cay constant of scalar meson 13 (1430), fx: FP?” (mke) 


is suppressed relative to fF? (m%). Therefore, the 


decay width of the first channel should be suppressed 
comparing with the second one. The numerical results of 
these two decay modes are given as below: 
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TABLE III: The branching ratios of B — K§(1430)K“) in the unite 1077. 


Decay Mode 
B” > Ko (1430)K° 
B- — K% (1430) K° 
B + K” K3°(1430) 
B- —> K` K#°(1430) 
B’ > K** (1430)K- 
DB’ + K+ KX (1430) 
B > Ko (1430) K*° 
B- — K% (1430) K*° 


DB? > K” KX? (1430) 


DB’ + K**(1430)K*~ 


11.121399- 5.57-5.05 
5.3611.37-2.36-2.30 
32.27 17°53! 6486 30 
23.7173 60 ae. 3.64 
+-0.43, 
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+0.07+0.07+0.09 
0.08" 9:02~0.04—0.06 
4 -0.21+0.05 
0.62 1008—0.18—0.04 
10.8673 34 
z EER +3.14+0.26+1.16 
B- + K* Kò (1430)| 12.717 575" 0.56" 1.07 
+0.14+0.00+0.22 
0.3310. 10—0.00—0.16 


0.084 


+0.25+0.53 
—0.24—0.75 


0.5810:29+0.03—0.05 
1.071047 -0.04-0.97 
0.030000.00—0.02 
0.14 75°05 "0.08" 0:03 
20.137 4.49"0.49"0.50 
21.7127 660.53-0.24 


+0.09+0.00+0.09 
0.111 0'06—0.00—0.06 


0.00+12.48 
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From these two equations, it is apparent that the emis- 
sion diagrams are dominant. However, there is a large 
enhancement from Og¢,g operators due to the fact that the 


chiral factor piso = 12.3 at u = 4.2 GeV is much larger 


than a = 1.5 owing to the larger mass of Kğ (1430). 


It follows that (ag — ža% yr is much greater than 
(ag — $ag)r and (aj — $a8)rX”. This compensate with 
the suppression of scalar meson decay constant to result 
in a larger branching ratio of B~ + K~ Kğ0 (1430) with 
scalar meson emitted. 

In the B => PP(V) and VV decay modes, the weak 
annihilation contribution is usually expected to be very 
small because it belongs to the next leading power correc- 
tion. However, one can see from Table III that the uncer- 
tainties induced by the weak annihilation are very large, 
even much larger than the central values in some decay 


modes, such as B° > Ky (1430). K 0 and pure annihila- 


tion mode B° > Ko (1430) K*—. This phenomenon can 
be understood as follow, when discussing the penguin- 
induced annihilation diagram of B + PP mode, the de- 
cay amplitude is helicity suppressed heavily because the 
helicity of one of the final-state mesons cannot match 
with that of its own quarks. On the contrary, this kind 
of helicity suppression can be alleviated when the scalar 
meson involved due to the nonzero orbital angular mo- 
mentum L, of the scalar meson. 

From Table III, one could notice that our predicted 
central values for the branching ratios of K?” (1430) K° 
and Kò (1430)K° based on S2 are smaller than the re- 
sults based on S1 by a factor of 4 ~ 5. In contrast, 


the central values of of K” K*(1430) and K~ K}°(1430) 


annihilation diagrams 


(32) 
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based on S1 are a bit smaller than those of S2. For 
K K0 (1430) and K*~ Kě0(1430), the predicted cen- 
tral values of S2 is larger than those of S1 by a factor 
2. Although there is difference between the central val- 
ues of different scenarios, we cannot distinguish two sce- 
narios due to large uncertainties taken by annihilation 
diagrams, unless one approach were proposed to deal 
with annihilations effectively in QCDF. In previous stud- 
ies [3], by comparing the calculated branching ratios of 
B — Kon and B > Kğ with experimental data, one 
concluded that the S2 is prefered, i.e. K is very likely 
the lowest lying state. If so, the branching fractions 
of K K%°(1430) and K~ K%(1430) could be measured 
by analyzing the data of K’ K+n- and K~K+a-. Un- 
fortunately, the predicted results of Ke (1430) K * and 
Ko (1430) K*° is too small to be measured with the cur- 
rent data of Belle, but they are hopeful to be measured 
in the forthcoming Belle-II by analyzing the four-body 
final states of KT Kat 17. 

Within the framework of PQCD, X. Liu et.al had cal- 
culated these decays [6], and both their results and ours 
are below the experimental upper limits. Comparing our 
predictions and theirs, we find that for pure annihilations 
they obtained rather larger branching fractions by keep- 
ing the transverse momenta, which are larger than ours 
by more than two order of magnitude. For decays with 
emission diagrams, they also got the larger branching ra- 
tios with rather large nonfactorizaion diagrams based on 
both two different scenarios. We hope the future mea- 
surements could distinguish two different frameworks. 

The CP asymmetries have not been observed in any B 
decays involving a scalar meson. For the charged decay 


TABLE IV: The direct CP asymmetries (%) of B > Kë (1430) K 


Decay Mode 
B- — Kj” (1430) K° 
B- > K~ Kö? (1430) 


+0.50+0.31+2.59 
—5.27 0.89—0.57—2.82 
=3 06+0-60+0.85+5.46 


—0.65—0.69—6.67 


S2 
+4.90+5.63+19.61 
— 22.5117 57—9.36—22.86 
+1.61+0.59+3.52 
= 2.60" 1 76~0.59—5.47 
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(33) 
the predictions of the direct CP asymmetries based on 
two different scenarios of scalar mesons are summarized 
in Table IV. The resource of this asymmetry is the inter- 
ference between the emission diagrams and annihilations. 
The large uncertainties in the latter lead to the large er- 
rors in the direct CP asymmetries, as shown in the table. 
For the neutral decay modes, the situation becomes 
more complicated due to the fact that both B° and B 
could decay to the same final states. For illustration, we 
will take B° (B°) 3 K?’ (1430) K° as an examples. Four 
decay amplitudes, Ay, A;, A f and A 7 are denoted as 


Ar = (K5K |B), 
= «770 50 
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A; =(Kj)K |B). (34) 
Then, we can write down the CP asymmetry as 
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Due to B° — B° mixing, the four time-dependent decay 
widths are given by (f = Ko’ (1430)K?) 


-rl z 
e ae Asl? + Ay 2) 
1 + Cy, cos Amt — Sp sin Amt], 


rE e> f) = ex (Az? +i 


(B?) + f) = 


2 
1- C cos Amt + Ssin Amt], 
1 = 
=P 
T(B%(t) + f) = e™ lA + Ap) 


1+ CF cos Amt — Sp sin Amt], 


— 1 = 
r(B (t) > f) = eTa(Ag?? +ã?) 
1 — Cy cos Amt + Sf sin Amt] (36) 


where Am stands for the mass difference of two mass 
eigenstate of B?°/B meson, and I for the average decay 
width of the B meson. The auxiliary parameters Cy and 


+1.30+0.17+11.40 
5.03 _1.37—0.17—13.52 


+2.54+0.23+7.57 
0.6413.64—0.22—9.22 


Sf appearing in above equations are defined by 


|As? — AFl? 
= E, (37) 
f JAHA 
2Im(A f) 
se = ——_~ 38 
t= DFAA (aa 
Va V Aş 
Ay = Dua, 39 
f= Wei Ay (39) 


Replacing f by f, we could obtain the formulaes for C F 
and Sp, correspondingly. If the experiment could find 
the values of Cf, Sp, Cy and SF, we thus can obtain four 
new parameters: 


1 
C = 5(C + On), AC = (Cs = CF), (40) 


wl ple 


1 
S= 3 (Sr + S7), AS =<=(Ss—S}). (41) 
Physically, S is the mixing-induced CP asymmetry and 
C is from the direct CP asymmetry, while AC and AS 
are CP-even under CP transformation Ay —> 1/A;7. In 
Table.V, we present our estimations of Acp, C, AC, S 
and AS for the final states Ky K°, Ko K”, KK- 
and KŠK *~ under two different scenarios. It should 
be stressed that in PQCD [6], there is no CP asymme- 
tries because of the absence of tree operators. However, 
in QCDF, by including the penguin contractions (P;), 
the charming penguin namely, an extra strong phase to- 
gether with another one from annihilations might lead to 
the large CP asymmetry, as shown in the table. Again, 
the large uncertainties of annihilations result in the large 
error for the pure annihilation decay modes, especially 
under S2. We hope these parameters can be measured in 
future colliders, such as high luminosity Belle-II, LHC-b 
and even higher energy ete” collider. 


IV. SUMMARY 


In this work, we studied the B > K (1430) kK“) decays 
under two different scalar meson scenarios by using the 
QCD factorization approach. We calculated the branch- 
ing fractions and the CP asymmetry parameters. It is 
found that the decay modes with the scalar meson emit- 
ted, have large branching fractions due to the enhance- 


* 


ment of large chiral factor re °, with some of the branch- 
ing fractions around the corner of Belle II. In contrast, 


TABLE V: The CP-violating parameters Acp, C, AC, S and AS (%) of B > K§(1430)K™. 
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S1 
S2 +0. 0040. “00-F0.00 


—*0 
Ko K° 


=*0 
Kọ? K“ 


S2 


+0.00-+0.00+0.00 
S1 —0.00—0.00—0.01 


—0.00—0.00—0.66 


the branching fractions of decay modes with the vector 
meson emitted, are much smaller. Moreover, the anni- 
hilation contributions take large uncertainties because of 
the free parameter from endpoint singularity. For the 
pure annihilation type decays, we predicted very small 
branching fractions, which are 2-3 orders of magnitudes 
smaller than the results from PQCD. Some of the decay 
channels are hopeful to be measured in future colliders, 
such as Belle-II, LHC-b and even high energy ete~ col- 
lider. 
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